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Abstract. —Despite its ancient origins, its worldwide distribution, and adaptation to diverse hab¬ 
itats, Isoetes has a highly conserved morphology. This feature has made it difficult to resolve 
species and species relationships using morphological characters. In this paper, we report the 
utility of nucleotide sequences from the nuclear internal transcribed spacer (ITS) regions, chlo¬ 
roplast atpB/rbcL intergenic spacer region, and second intron of a LEAFY (LFY1 homolog for iden- 
tdying species relationships, delimiting basic diploid species, and determining hybrid origins. 
Variation in the ITS regions and atpB/rbcL spacer is most useful at the family level in Isoetes and 
the LF 1 second intron is appropriate at the species and population level. The tree resulting from 
an analysis of the combined nuclear ITS and chloroplast ntpB/rbcL spacer contains three major 
well supported clades (bootstrap > 99%): an Old-World/California clade (/. abyssinica, I. longis- 
sima. I. velata, I. nuttallii, and I. orcuttii), an Asian/Australian clade (/. taiwanensis, /. japonicci, 
I. kirkn. and 1. drummondii), and a poorly resolved clade consisting of nine North American 
species. To further resolve and delimit the North American species, a combination of the LEAFY 


second intron and ITS data was used. The resulting consensus tree has limited resolution, sup¬ 
porting the hypothesis that the North American species complex radiated rapidly. The combina¬ 
tion of LEY and ITS data provided numerous characters, both substitutions and indels, that are 


useful in species delimitation and identification of cryptic species. ITS sequence data, through 

additive banding and sequence misalignment, is also useful in confirming interspecific hybrids 
and determining their parental origins. 


Isoetes is a cosmopolitan genus of heterosporous lycopsids containing 200 
or more species. Despite the probable ancient Paleozoic origins of Isoetes. its 
worldwide distribution, and its diverse ecological adaptations (from seasonal 
ephemeral terrestrials to obligate evergreen aquatics), Isoetes species are sin¬ 
gularly lacking in morphological variation (Taylor and Hickey 1992). The tax¬ 
onomic difficuties caused by conserved morphology are compounded by ho- 
moplasy and allopolyploid speciation. Because of these problems. Isoetes is a 
prime candidate for additional sources of data, such as DNA sequences. 

A typical Isoetes sporophvte consists of a buried rootstock bearing a rosette 
of simple, linear sporophylls. Each sporophyll contains four longitudinal, sep¬ 
tate canals and an ovoid sporangium near its base. Each sporangium contains 
either scores of globose megaspores or thousands of ellipsoid microspores. A 

flap ol tissue, called the velum, partly or completely covers the sporangium. 
Above the sporangium is a triangular ligule. 

Since the Paleozoic Era, there have been numerous opportunities tor species 
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of Isoetes to adapt repeatedly to terrestrial habitats during drier conditions and 
to aquatic habitats in wetter periods. The range of habitats in which the species 
are found today appears to be a reflection of these past adaptations. For ex¬ 
ample, some species are evergreen aquatics in glacial lakes. Other species are 
ephemeral, seasonal terrestrials in soil pockets on exposed rock outcrops. Still 
others are amphibious, found in and along intermittent streams. A few species 
invade old fields and roadside ditches. Isoetes appears to be an opportunistic 
pioneer able to adapt to a wide range of habitats. 

Repeated adaptations to aquatic and terrestrial environments over time 
could lead to reversals, parallelisms, and convergences in morphological data 
sets. Thus it is possible that homoplasious morphological character states have 
been used to identify and relate species. A phvlogeny of Isoetes, based on lew, 
potentially homoplasious character states, would probably be an inaccurate 
reconstruction. Therefore, alternative data sets are critical for reconstructing 
phylogenies and defining species within Isoetes. 

Chromosome counts have provided some insight into the evolutionary his¬ 
tory of Isoetes. For example, about one-half of the North American species are 
basic diploids (2n = 22, X = 11). The other half form a polyploid series of 
tetraploid (2n = 44), hexaploid (2n = 66), octaploid (2n = 88), and decaploid 
(2n = 110) species. These chromosome numbers indicate that Isoetes has 
evolved in two ways. Basic diploid species have evolved by ecological isola¬ 
tion and genetic divergence as taxa adapted to terrestrial and aquatic habitats. 
Polyploid species have evolved through interspecific hybridization and chro¬ 
mosome doubling when divergent species occupied overlapping ranges. Inter¬ 
specific hybrids, recognized by their irregularly formed and nonviable mega- 
spores, are commonly encountered where two or more species occur together. 
At least fifteen different interspecific hybrids of Isoetes have been detected 
among North American species (Taylor et al. f 1993; Montgomery and Taylor, 
1994; Britton and Brunton, 1995, 1996; Musselman et al., 1995, 1996, 1997). 


Interspecific hybrids and putative allopolyploids have caused confusion in 
defining species because they possess character states intermediate between 

the two parental species. 

The main purpose of this paper is to show the usefulness of sequence data 
for ITS. the citpB/rbcL spacer region, and especially the second intron of the 
LEAFY ( LFY) ortholog, to address phylogenetic questions. Future work will 
result in more extensive data (including many more species sampled) for re¬ 
solving phylogenetic relationships and questions related to biogeography, spe¬ 
cies delimitation, and hybrid origins. 

Variation in the ITS regions and citpB/rbcL spacer is most useful at the family 
level in Isoetes and the LFY second intron is appropriate at the species and 
population level. The two ITS regions are located between the 18S and 26S 
nuclear ribosomal genes and separated b\’ the 5.8S gene. There are multiple 
copies of this ribosomal array' in the genome, but they appear to undergo rapid 
concerted evolution and all copies appear to be virtually identical in Isoetes 
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insertions/deletions (indels) are common. ITS regions have been used exten¬ 
sively to study angiosperm taxa at the species and generic level (Baldwin et 
al, 1995). 

The atpB/rbcL intergenic spacer is located between two highly conserved 
chloroplast genes, atpB and rbcL. This spacer region is approximately 750 bp 
long, but indels are common. Spacer sequences have been used in many phy¬ 
logenetic studies (e.g., Golenberg et al., 1993; Manen et al., 1994; Hoot and 
Douglas, 1998). 

LFY is a meristem identity gene that encodes a transcription factor involved 
in the activation of homeotic genes which control flowering in Arabidopsis 
(Parcy et al., 1998). LFY homologs have been identified in Pinus, Picea, Wel- 
witschia, Gnetum, ferns, Huperzia, Lycopodium, and Selaginella (Frohlich and 
Meverowitz, 1997; Mouradov et al., 1998; Frohlich, pers. comm.; Therrien, 
pers. comm.). LFY is usually found in one copy in higher plants but Isoetes 
has two copies. These two copies, including short stretches of the flanking 
exons, differ in length by about 100 bp (long copy ca. 1150 bp; short copy ca. 
1050 bp). The partial exon sequences of the two copies are highly conserved 
and readily aligned, but the two introns have very different nucleotide se¬ 
quences. For our work on Isoetes, only the long copy of intron 2 was se¬ 
quenced. 


Materials and Methods 

Total DNA was isolated from 0.3—2.0 g of fresh or silica-dried leaf tissue 
(Table 1) following the CTAB extraction procedure of Doyle & Doyle (1987). 

The ITS regions, including the 5.8S gene, were amplified by the polymerase 
chain reaction using primers 1830F, located in the 18S gene, and 25R, located 
in the 26S gene (for primer sequences, see Nickrent et al., 1994). Double- 
stranded amplifications in 100 jjlL reactions were conducted with the following 
reagents: 20 raM Tris (pH 8.3), 50 mM KC1, 1.5 mM MgCl.,10 M DMSO. 50 
jjlM each dNTP, 0.25 jxM of each amplification primer, and 2.5 U Taq poly¬ 
merase. After overlaying each reaction with mineral oil, genomic DNA (210- 
500 ng) was added. Reaction conditions consisted of 40 cycles of 94° C for 30 
sec, 45°C for 30 sec, and 72 C for 2 min for the denaturation, annealing, and 
extension steps, respectively. The first cycle was preceeded by a 4 min dena¬ 
turation step and the last cycle by a 5 min extension step. The atpB-rbcL in¬ 
tergenic spacer region was amplified as described in Hoot and Douglas (1995) 

using primers located within the atpB and rbcL genes for both amplification 
and sequencing. 

The second intron of LFY (long copy) was amplified using the degenerate 
primers LFsxl-3 and LFtxr (Frohlich and Meyerowitz, 1997). For five species, 
the two copies of LFY present in Isoetes were substantially different in size 
and could be separated on a 1.5% agarose gel overnight. The two bands were 
then excised, reamplified using the melted gel plug as template, and used for 

sequencing. The resulting sequences were then aligned to locate a conserved 
region within the second intron that could ha used in datiion intomol amrUi. 
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Table I. Collection data tor species analyzed. 


Isoetes 

Country/State 

C ollector 

Collection 

/. abxssinica Chiov. 

Kenya 

Gastony 

97-101 

/. bu fieri Engelm. 

USA/Kansas 

McGregor 

s.n. 

/. coromandelina L. 

India 

Srivastave 

SM. 

/. d rum monel ii A. Br. 

Australia 

Hoot 

s. n. 

/. echinospora Dur. 

USA/Maine 

Taylor 

5137 

/. echinospora Dur. 

USA/Montana 

Taylor 

s.n. 

/. echinospora Dur. 

USA/Wisconsin 

Taylor 

5097 

/. engelmannii A. Br. 

USA/Michigan 

Taylor 

s.n. 

/. Engelmannii A. Br. 

USA/Missouri 

Taylor 

s.n. 

/. Engelmannii A. Br. 

USA/Virginia 

Taylor 

6015 

/. flaccida Shuttlew. 

USA/Florida 

Taylor 

5214 

/. hawaiiensis Taylor & Wagner 

USA/Hawaii 

Taylor 

5658 

/. japonic a A. Br. 

Japan 

Takahashi 

s.n. 

/. kirkii A. Br. 

New Zealand 

Woodland 

s.n. 

/. lithophila Pfeiffer 

tJS A/Tex as 

Taylor 

4653 

/. longissima Bory 

Spain 

Taylor 

5409 

/. melanopoda Gay & Dur. 

USA/Louisiana 

Lark 

s.n. 

/. melanopoda Gay & Dur. 

USA/Mississippi 

Leonard 

s.n. 

/. melanopoda Gay &. Dur. 

USA/Arkansas 

Johnson 

s.n. 

/. melanospora Engelm. 

USA/Georgia 

Taylor 

4849 

/. nuttallii A. Br. 

USA/Califomia 

Taylor 

5446 

/. orcuttii A* A. Eat. 

USA/Califomia 

Taylor 

5447 

/. prototypus Britton 

Canada 

Brunton 

12350 

/. taiwanensis DeVol 

Taiwan 

Kuo 

s. n. 

/. te get iformans Rury 

USA/Georgia 

Taylor 

5182 

/. valida (A. A, Eat.) Clute 

USA/North Carolina 

Taylor 

s.n. 

/. velata A. Br. 

Spain 

Prada 

s.n. 


fication primers. These internal primers were specific to just the long copy of 
the LFY intron and are located approximately half way through the intron. 
The long copy of the LFY second intron was then amplified in two parts. 
Amplification followed the same protocol described above tor ITS except that 
an initial 25 p.1 amplification reaction was done with Ready-To-Go Beads 
(Amersham Pharmacia Biotech) and the resulting PCR product used as the 
template for future reactions. 

Following amplification, double-stranded PCR product was treed ot remain¬ 
ing nucleotides and primers by electrophoresis on a low melting 2% agarose 
gel. The bands were visualized with UV light source, removed as gel plugs, 
and melted at 70°C. DNA was further purified and concentrated using Wizard 
Columns (Promega). following the manufacturers protocols. Purified PCR 
products were sequenced on an ABI 373-Stretch automated sequencer. The 
two amplification primers were used tor each region with additional internal 
sequencing primers employed for the ITS and atpB-rbcL spacer regions (primer 

sequences available upon request from S. Hoot). 

Alignment of DNA sequences was accomplished to a rough approximation 
using Sequencher 3.0 (Gene Codes Corp.) with subsequent manual corrections. 
The criteria for indel alignment and scoring are as described in Hoot and Doug- 
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las (1998). Phylogeneticallv informative indels (variable in two or more taxa) 
were scored as one event at the end of the data set. 

Fitch parsimony of the separate and combined data sets were performed 
with PAUP* vers.4.0b2a (Swofford, 1998) using the branch and bound search 
option. To assess branch support. PAUP* was used to perform bootstrap anal¬ 
yses (Felsenstein, 1985) with 500 replications, the branch and bound search 
option, and maximum number of trees saved set at 2000. 

Isoetes is morphologically very distinct and does not appear to be closely 
related to any extant pteridophytes. We are currently working on rooting our 
molecular trees using various genera within the Lycopodiaceae, the nearest 
sister group to Isoetes (Manhart, 1994; Wikstrom and Kenrick 1997). As this 
work is not yet complete, we chose I. coromandelina L. f. as the root of our 
tree based on three lines of evidence: 1) preliminary rbcL data for selected 
species of Isoetes and lycopods, 2) previous morphological work suggesting 
that the more piesiomorphic character states and geographic distributions are 
found in section Coromandelina (Taylor & Hickey, 1992), and 3) molecular 
evolution in Isoetes appears to evolve in a clocklike fashion, suggesting that 
the root of the tree will be among the species with the greatest number of 
substitutions (this work and unpublished preliminary data). For the analysis 
of the combined LFY/ITS data (North American species complex), the tree was 
rooted using the ITS sequence for I. drummondii A. Br., based on the results 
of our worldwide phylogeny (Fig. 1). Due to sequencing difficulties, a sequence 
for I. drummondii was not available for LFY ; the data for that region is scored 
as missing for this species in the combined data set. 

Results 

Relationships of Basic Diploid Species. —A combination of ITS and atpB-rbcL 

spacer data was used to evaluate relationships on a global scale and resulted 
in seven most parsimonious trees based on 176 informative characters (312 
variable characters) with a consistency index excluding uninformative char¬ 
acters (Cl) = 0.79 and retention index (Rl) = 0.90. One of the seven shortest 
trees is presented in Fig. 1. 

Three major well supported clades (bootstrap ^ 99%) are recognized: an 
Old-World/California clade (/. nuttalii, /. orcuttii, /. abvssinica. /. longissima. 
and /. velata), an Asian clade (/. kirkii, I.drummondii, /. taiwanensis. and /. 
japonica), and a North American clade (/. flaccida, I. prototypus. /. melano- 
poda, I. lithophila, /. melanospora, /. engelmannii, I. valida. I. echinospora, 
and /. hawaiiensis ). The North American species are not monophyletic. Isoetes 
nuttallii and /. orcuttii from California are nested within a clade consisting of 
Spanish and African species, separate from the remaining North American 
species. The North American clade, referred to in the future as the North 
American species complex, is poorly resolved, suggesting a possible rapid ra¬ 
diation event with subsequent speciation. 

Species relationships and delimitation within the North American spe¬ 
cies complex. —To resolve and delimit species within the North American spe- 
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coromandelina India 
nuttallii CA 
orcuttii CA 
abyssinica Kenya 
longissima Spain 
velata Spain 
kirkii New Zealand 
drummondii Australia 
taiwanensis Taiwan 
japonica Japan 
valida VA 

echinospora ME/WI/MT 
hawaiiensis HA 
melanospora GA 
engelmannii MI/MOA/A 
flaccid a FL 
prototypus NB 
melanopoda LO/MS 
lithophila TX 


Fig. 1 . One of the seven most parsimonious trees derived from the combined ITS and atpB-rbcL 
spacer region data. Numerals above the lines indicate the number of nucleotide changes support¬ 
ing each branch (using accelerated transformation optimization); numbers below branches are 
bootstrap values. Dotted lines indicate branches that collapse in the strict consensus tree derived 

from multiple shortest trees. 


cies complex, a combination of ITS and LFY intron data was used. The LFY 
intron data, consisting of 218 variable characters and 24 indels, is more than 
four times more variable than ITS (48 variable characters and four indels). The 
eight trees resulting from the analysis of LF\ intron data alone (tree not pre¬ 
sented) are largely similar to the combined LFT /ITS tree (Fig. 2), differing onl\ 
in the position of /. butleri. (sister to /. engelmannii in the LFY tree). There 
was no resolution of the North American species complex using ITS data 

alone. 
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engelmannii MI 


engelmarmii MO 


engelmannii VA 



melanopoda MS 


melanopoda LA 



hawaiiensis HA 


echinospora MT 


echinospora ME 


echinospora WI 

Fie. 2. One of six most parsimonious trees derived from the combined LFY intron 2 and ITS 

d<it«i. Numerals above the lines indicate the number oi nucleotide changes supporting each branch 

(using accelerated transformation optimization); numbers below branches are bootstrap values. 

Dotted lines indicate branches that collapse in the strict consensus tree derived from multiple 
shortest trees. 


One of the seven shortest trees (94 informative characters. Cl = 0.77 and RI 

= 0.88) resulting from the analysis of the combined ITS and LFY intron data 

is presented in Fig. 2. The deeper branches of the tree collapse in the strict 

consensus tree, supporting the hypotheses that the North American species 
complex radiated rapidly. 

Within the North American species complex, there is some resolution of 
species relationships (Fig. 2). Isoetes tegetiformans, I. lithophila. and I. melan¬ 
opoda (population from Arkansas) form a moderately well supported clade (9 
substitutions, bootstrap = 71%). Unlike the results from a combination of ITS 
and spacer data, I. melanopoda populations from Louisiana and Mississippi 
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are weakly supported (bootstrap = 63%) as sister to a well supported clade 
consisting of 7. hawaiiensis and 7. echinospora (bootstrap = 100%). While the 
position of 7. valida changes between the two trees derived from ITS/spacer 
and LFY/ ITS, none of the relevant branches are well-supported (bootstrap < 
58%). 

The combination of LFY' and ITS data provide numerous characters, both 
substitutions and indels, which are useful in delimiting species (Fig. 2). Plants 
from three widely separated populations of 7. engehnannii from Michigan, Mis¬ 
souri, and Virginia occur in the same well-supported clade (Fig. 2; 32 char¬ 
acters, bootstrap = 100%). Similarly, individuals from widely separated pop¬ 
ulations of I. echinospora from Maine, Wisconsin, and Montana, are nearly 
identical, only differing from each other by three unique sites found in the 
Montana population. 

In contrast, the tree resulting from the combined LFY' and ITS data indicates 
that 7. melanopoda as currently defined is polyphyletic and composed of at 
least two species. Isoetes melanopoda from Arkansas, differs from two popu¬ 
lations of 7. melanopoda from Louisiana and Mississippi by numerous substi¬ 
tution and indels. The former appears in a moderately-supported clade with 
I. lithophila and /. tegetiformans (bootstrap = 71%). 

Determining the Origin of Allopolyploid Species. —We used nucleotide se¬ 
quences of the ITS regions of nuclear ribosomal DNA to tentatively determine 
the origins of the interspecific hybrid /. Xeatonii Dodge and the allotetraploid 
7. riparia Engelm. ex A. Br. Taylor and Hickey (1992) provided data from dis¬ 
tribution patterns, spore morphology and viability, electrophoretic profiles of 
leaf enzymes, and in vitro hybridization experiments supporting the hypoth¬ 
esis that the basic diploid species 7. echinospora (2n = 22) and 7. engehnannii 
(2n = 22) have crossed to form the sterile hybrid 7. Xeatonii (2n = 22) and 
that chromosome doubling has produced the fertile allotetraploid, 7. riparia 
(2n = 44). 

When ITS sequences of 7. echinospora, I. engehnannii, their hybrid, 7. xea¬ 
tonii, and 7. riparia are aligned and compared, nucleotide additivity and mis¬ 
alignment are apparent in the sequence chromatograms (Fig. 3). Nucleotide 
additivity occurs in 7. xeatonii and 7. riparia at nucleotide positions where 7. 
echinospora and 7. engelmannii differ. Nucleotide misalignment occurs in 7. 
Xeatonii and 7. riparia in the 5' to 3' direction (determined by the direction 
of primer extension) where there is a nucleotide insertion or deletion in either 
7. echinospora or 7. engelmannii (Fig. 3). ITS nucleotide additivity occurs at 
each of 10 sites in 7. Xeatonii and at 8 sites in 7. riparia where 7. echinospora 
and 7. engehnannii differ (Table 2). There are also two indels in 7. echinospora 
and engelmannii, leading to nucleotide misalignments (beginning at sites 49 
and 548) in 7. Xeatonii and 7. riparia. 

Discussion 

Relationships of Basic Diploid Species.— The combination of nucleotide se¬ 
quences from ITS and the atpB/rbcL intergenic spacer are effective in resolving 
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i<- 

Cfi RTC ft C 6 7 G G T C 6 6CTG 



Cfl flGC flCITC G TCG GC1G 


Isoetes X eatonii 



Cfl RttCflCT HGTCHGHT6 



Cfl R8C flCHTGHHCe GHTG 


Fig. 3. Aligned sequence chromatograms of /. engehmmnii , /. echinosporci, /. x eatonii. and /. 
viparia for an 18 bp segment in the 1TS2 region. The arrow on the left indicates site 544 where 
the two parental species, /. engehnannii and /. echinospora . differ in nucleotides and additive 
banding is found in the two taxa of hybrid origin, /. Xeatonii. and /. riparia. The arrow on the 
right indicates the location of a deletion in /. echinospora in relation to /. engehnannii. The two 
parental sequences found in /. Xeatonii and /. riparia are misaligned to the right of this deletion 
in the direction of the primer extension. 


phylogenetic relationships of diploid species found throughout the world. 
These variable sequences and their rate of evolutionary change provide nu¬ 
merous characters (both nucleotide substitutions and indels) at this taxonomic 
level in Isoetes. 

The tree resulting from the combination of ITS and spacer data are of interest 
biogeographically (Fig. 1). The well supported Old World/California clade con¬ 
sists of species with Eurasian affinities. The inclusion of the two Californian 
species. I. nutallii and /. orcuttii, within this clade suggests that this group of 
Isoetes may have previously had a distribution throughout Asia and into North 
America. Isoetes abyssinica is sister to the two Spanish species, /. longissima 
and /. velata, indicating a strong affinity between at least some North African 
and European species. These results also support at least two origins for the 
North American species of Isoetes, those with Eurasian affinities as described 
above and the many species found in the North American species complex. 
Current work is now underway which will greatly increase our sampling 
worldwide, further testing the above biogeographical hypotheses. 

The lack of variation in the North American species complex using two 
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Table 2. ITS nucleotide additivity in Isoetes X eatonii and /. riparia at sites where the parental 
species, /. englemannii and /. echinospora , differ. 


Taxon 




Sites and Nucleotides 




39 

83 

274 

314 

319 

321 

544 

627 

644 

679 

I. engelmannii 

A 

C 

T 

G 

G 

A 

T 

c 

G 

G 

/. echinospora 

T 

T 

c 

C 

T 

T 

G 

T 

A 

T 

/. X eatonii 

Art 

C/T 

C/T 

C/G 

G/T 

ATT 

G/T 

C/T 

AJG 

G/T 

L riparia 

T 

Crt 

C 

C/G 

G/T 

Art 

G/T 

Crt 

A/G 

G/T 


different combined data sets, ITS/spacer (Fig. 1) and LFY7ITS (Fig. 2), suggests 
a rapid radiation in the North American species complex. However, the LFY/ 
ITS tree does resolve a sister group relationship between /. hawaiiensis and 
three populations of /. echinospora (8 substitutions, 2 indels, bootstrap = 
100%), suggesting that a founder population of I. echinospora made its way 
to Hawaii from the boreal regions of the Northern Hemisphere, possibly by 
means of migrating waterfowl. Isoetes hawaiiensis, while exhibiting consid¬ 
erable morphological divergence from /. echinospora, is virtually identical in 
terms of molecular data (two character differences), suggesting that the event 
occurred relatively recently. 

•/ mJ 

Species delimitation within the North American species complex. —A per¬ 
sistent problem in molecular plant systematics has been the lack of DNA se¬ 
quence regions variable enough to be used at the species and population level. 
The second LFY intron is proving invaluable in this regard. In Isoetes, LFY is 
at least four times more variable than ITS, which is generally regarded as one 
of the more variable sequence regions available for systematic studies. 

The species identity found in the LFY/ ITS tree (Fig.2) is largely based on 
LFY data due to the greater variation found in this region. LFY provides nu¬ 
merous characters, both substitutions and indels. which can serve to delimit 
and identify species. For example, our results indicate that widely separated 
populations of I. engelmannii (Michigan, Missouri, and Virginia) and I. echin¬ 
ospora (Maine, Montana, and Wisconsin) are virtually identical. 

In contrast, the I. melanopoda complex contains at least one cryptic species. 
The population of /. melanopoda from Arkansas is very different molecularly 
from the two populations from Louisiana and Mississippi, coming out within 
a moderately supported clade (bootstrap = 71%) consisting of I. lithophila 
from Texas and /. tegetiformans from Georgia. Work is currently underway to 
further test species limits for other Isoetes taxa. 

Determining the Origin of Allopolyploid Species. —Additive banding and 

sequence misalignment in ITS data from hybrids and allopolyploids are valu¬ 
able tools in the confirmation of interspecific hybrids and the determination 
of their parental origins. For the most part, the additive banding found in the 
hybrid /. X eatonii and the allotetraploid /. riparia confirm that their progeni¬ 
tors were /. echinospora and /. engelmannii (Table 2). Two of the sites (39, 
274) do not show additivity in /. riparia, but do possess nucleotides at these 
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sites that are found in one of the putative parents. This loss of additivity may 
be due to the effects of concerted evolution, mutation events, or sequence 
differences in one of the diploid parents of I. riparia not detected with our 
sampling. We are currently studying hybrid origins of the above and other 
allopolyploids using LFY data. Because many more substitutions and indels 
occur with LFY data, cloning is being employed to separate the two parental 
copies. Our preliminary data for LFY also support the above findings for the 
origins of I. riparia. 
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